Solution of the Stochastic Langevin Equations for Clustering of 
Particles in Random Flows in Terms of Wiener Path Integral 

M. Chaichiani'^y A. Tureanu^-^lll and A. Zahabii 

^Department of Physics, University of Helsinki, 
P.O. Box 64, FIN-OOOI4 Helsinki, Finland and 
^ Helsinki Institute of Physics, P.O. Box 64, FIN-OOOI4 Helsinki, Finland 

(Dated: February 7, 2010) 

Abstract 

We propose to take advantage of using the Wiener path integrals as the formal solution for the 
joint probability densities of coupled Langevin equations describing particles suspended in a fluid 
under the effect of viscous and random forces. Our obtained formal solution, giving the expression 
for the Lyapunov exponent, i) will provide the description of all the features and the behaviour of 
such a system, e.g. the aggregation phenomenon recently studied in the literature using appropriate 
approximations, ii) can be used to determine the occurrence and the nature of the aggregation 
- non- aggregation phase transition which we have shown for the one-dimensional case and Hi) 
allows the use of a variety of approximative methods appropriate for the physical conditions of 
the problem such as instanton solutions in the WKB approximation in the aggregation phase for 
the one-dimensional case as presented in this paper. The use of instanton approximation gives 
the same result for the Lyapunov exponent in the aggregation phase, previously obtained by other 
authors using a different approximative method. The case of non-aggregation is also considered in 
a certain approximation using the general path integral expression for the one-dimensional case. 
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I. INTRODUCTION 



Recently, there has been considerable interest in the study of the behaviour of particles in 
media taking into account the effect of random forces. Studies in this direction can provide 
a better understanding of the behaviour of particles in turbulent flow. The features and 
behaviour of turbulent flow are under continuous intense investigations. As one of the main 
features, the clustering of particles into regions of high density has been studied extensively 
on both experimental and theoretical sides . Particles suspended in a turbulent fluid 
form cluster structures as a result of the competition between the diffusive random forces 
and the aggregative viscous ones. However, the conditions for such a behaviour are not fully 
understood, the mechanisms which contribute to the formation of clusters have been studied 
m (see also j9|). An extreme form of clustering of particles, known as the "aggrega- 

tion phenomenon", which is not well-understood, has been studied recently by means of 
theoretical modeling and numerical simulations 10|]- 16|. Other phenomenological models 



for cluster aggregation, inspired by Kolmogorov's theory 17|], may be also studied along 
similar lines [la]. The aggregation of particles can be defined as the coalescence of different 
particles paths with very close positions and velocities in a fluid subjected to random forces 
fluctuating in space and time, the particles being affected by viscous forces proportional to 
their velocities. The first theoretical analysis and numerical simulations for th e ag gregation 
of suspended particles in a one- dimensional random fluid were carried out in 10|. The re- 
sult of this study shows a phase transition between the non-aggregate and aggregate phases. 
Motivated by this result, recent investigations on the ag greg ation of particles in two- and 
three-dimensional random fluids were performed in jll|-|l3|. By introducing a model for 
the motion of point-like non-interacting particles in a three dimensional random fluid, the 



equations of motion for such particles, w. 
the random force were also derived in 



lich are under the influence of a viscous force beside 



13| . Then by linearizing the equations of motion. 



two coupled Langevin equations which describe the evolution of the separation of positions 
and velocities of two nearby particles were obtained. These two coupled Langevin equations 
describe the aggregation of particles and thus the system of coupled Langevin equations 
should be solved for calculating the Lyapunov exponent [19|, which is equal to the expecta- 
tion value of one of the variables in the Langevin equations. For this purpose, we apply the 
Wiener path integral formalism for solving the system of two coupled Langevin equations. 
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describing the aggregation phenomenon. At first we introduce a method for writing the solu- 
tion of Langevin equations in terms of the Wiener path integral, which has been thoroughly 
studied in the literature over the years 20|- 2J]. Then, by generalizing the procedure, we 
obtain the solution of the system of N coupled Langevin equations in terms of the Wiener 
path integral. The Lyapunov exponent as an indicator of aggregation can also be written 
in terms of path integral. The Wiener path integral formalism provides an exact solution 
to the aggregation problem. The obtained exact solution in terms of the path integral is 
presented in a closed analytical form and its actual evaluation can be performed by means 
of a variety of different approximation methods suitable to the specific physical conditions 
of the system. We discuss some approximative methods for such Wiener path integrals in 
one-dimensional case. 



II. SYSTEM OF COUPLED LANGEVIN EQUATIONS DESCRIBING THE AG- 
GREGATION PHENOMENON 



In this section we briefly review the results of recent studies [13| on the aggregation 
phenomenon in three dimensions. Particles suspended in a turbulent fluid can be modeled 
by the spherically massive particles which are moving in a random velocity fleld with speciflc 
properties such as isotropic, homogeneous, and stationary statistics. For simplicity, it can 
be assumed that there is no interaction between the particles themselves as well as between 
the particles and the fluid. Also we can neglect the inertia of the displaced fluid. By these 
assumptions, one can consider a large number of suspended particles with random initial 
positions in the fluid and zero velocities. The behaviour of the trajectories of particles is 
dictated by the effect of random as well as viscous forces and the motion of such particles in 
a random fluid is diffusive. Therefore, the inhomogeneities in density tend to get reduced, 
while the viscous forces cause the aggregation of particles and eventually the competition 
between diffusive random forces and viscous forces leads to a phase transition between path 
coalescence and path non-coalescence phases. The equations of motion which describe the 
suspended particles' trajectories are 

r = — , p = -7[p-mu(r,t)], (1) 
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where 7 characterizes the strength of the viscous damping and u(r, t) is the random velocity 
field. The aggregation of particles can be studied by considering two nearby trajectories with 
spatial separation Sr and momenta difference 5p. The linearized version of the equation of 
motion can be derived as 

6i = ^, 5p = -75p + F(t)5r, (2) 

m 

with the matrix elements of F as 

F., = §-m,t)=rnj^^{r{t),t). (3) 

With the parametrization of the linearized equations of motion as 

6r = Xm, 6p = X(rini + Fans), (4) 

where X is a scale factor showing the aggregation occurs if it decrease with probability in 
the limit t 00. ni and 112 are time-dependent orthogonal unit vectors (in three dimensions 
the third unit vector is defined as ns = rii A n2). By differentiating eq. (4), substituting the 
result in eq. (2) and finally projecting the equations onto the unit vectors one obtains the 
equations of motion for the variables X and 

X = —YiX 
m 

m 

% = -jY2--Y,Y2 + F^,{t), (5) 
m 



where nj(t).F(t)nj(t) = F'-{t). It has been argued [13|, that the maximal Lyapunov expo- 
nent Ai (hereafter called simply Lyapunov exponent) in the limit t — ?■ 00 is given by 

,d\iiX , 1 , , 

which we used the first equation in (5). 

We recall the general definition of the Lyapunov exponent as a quantity that characterizes 
the rate of separation of infinitesimally close trajectories. Quantitatively, the separation 
5Z{t) of two trajectories in phase space with initial separation 5Zq is given by the formula 

I 5Z(t) 1^1 5Zo I exp(At), (7) 
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where A is the Lyapunov exponent. The separation diverges with time when A > and 
aggregation occurs for the case A < 0. 

In the hmit where the correlation time r' of the random force is small and the random 
force itself is also sufficiently weak, the coupled equations of motion can be approximated 
by a system of two coupled Langevin equations (for details of derivations and physical 
explanations see part III of [13]): 



dYi 
dY2 



-^Y, + -{Yi-Y,' 
m 

-7^2 + -T7 Y1Y2 

Y2 m 



dt + dC^i, 
dt + dC,2-, 



with the noise properties 



(8) 



{dQ = 0, 



{dQdQ = 2D,jdt. 



(9) 



The diffusion constants Dij is defined as 



Consequently, by the change of variables 



dt{F:,it)F'iO))- 



dt' = 'fdt, 



Xi 




—Y- 



where Di = Du, D2 = -D21 = -D31, the final coupled Langevin equations read as 

dxi = [—xi + e{Txl — xf)] dt' + dwi, 
dx2 = [—X2 + X2'^ — 2^X1X2] dt' + dW2, 



(10) 



(11) 



(12) 



where e = /m'-f^^'^ is a dimensionless measure of the inertia of the particles and F = 
D2/D1 is a measure of the relative intensities of potential and solenoidal components of the 
velocity field. Equivalently, the coupled Langevin equations can be written as 



xi = [-xi + e{Txl - xj)] + liii, 
X2 = [-X2 + X2 ^ - 26X1X2] + W2, 



(13) 



with 



(wi) = 0, 



1,2. 



(14) 
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This system of stochastic differential equations should be solved and the derived solution, 
which represents the probability density, can be used to determine the Lyapunov exponent 
Ai [l9|, [l3|, as an indicator of the aggregation phenomenon: 



Xi = ie{xi). (15) 

When the Lyapunov exponent Ai is negative the aggregation phenomenon occurs flS^. A 
positive Lyapunov exponent is an indication that the system is chaotic. 

The system of two coupled Langevin equations as a special kind of stochastic differential 
equations has an exact solution in terms of the Wiener path integral. In the following, we 
introduce the Wiener path integral as the solution of stochastic differential equations for 
the special case of the Langevin equations. Subsequently, we generalize the Wiener path 
integral formalism for the system of coupled Langevin equations, which describe massless as 
well as massive Brownian particles in random media. After that we shall be able to present a 
solution in terms of the Wiener path integral for the coupled Langevin equations describing 
the aggregation phenomenon. 



III. SOLUTION OF THE STOCHASTIC DIFFERENTIAL LANGEVIN EQUA- 
TIONS IN TERMS OF WIENER PATH INTEGRAL 

In this section, we introduce the Wiener path integral method for the solution of stochastic 
differential equations of special kind, the Langevin equations. As a prototype for such 
Langevin equations, we can look upon them as describing the Brownian motion in different, 
coordinate or velocity, spaces. For a comprehensive description of several Brownian particles 



in general, see 



25|. 



A. Wiener path integral for one Langevin equation 

At first, we consider a Brownian particle in a random medium and write the path integral 
solution for the transition probability of the particle from one fixed arbitrary initial point to 
a fixed arbitrary final point. The Wiener path integral method can be used for the analysis 
of the stochastic equations and consists in determining the statistical properties of their 
solutions such as probability densities and expectation values. The microscopic approach to 
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stochastic processes starts from the stochastic Langevin equation. The Langevin equation 
for a Brownian particle subject to a general non-stationary and nonlinear external force is 



mx -\- r]x — F -\- ^ , 



(16) 



where m is the mass of the particle, r] is the friction coefficient, F is an external force and 
$ is a random force. For sufficiently large time intervals t m/r] we can neglect the mass 
term, so the Langevin equation which describes the motion of inertialess Brownian particles 
takes the form 

x{T)+f{xiT),T)=<P{T), (17) 

where 

F ■ $ 

/=— , </> = -. (18) 

—T) T] 

Performing a functional change of variables through the Volterra integral equation 

y{r) = x{T)+ ! f{x{s),s)ds, 0<T<t, (19) 
Jo 

one can write the eq. (17) as 

^(r)=0(r). (20) 
The Jacobian of this transformation can be evaluated by the discrete-time approximation 



l + f{x,,e)l 

f'{xi,s) 1 + /'(X2,2£)| 








f{x^,e) f{x2,2e) l + f{xr,,Ne) 

where /' = df /dx and e — t/N. The determinant (21) becomes 

\ ^ ef'{xn,ne) 



(21) 



N 



Jie) = n 



n=l 



(22) 



In the continuum limit, the determinant takes the form 



J = limexp 



£->0 



exp 



1 ^ 



n=l 



\jj\x{s),s)ds 



(23) 



Now we can write the transition probability of the stochastic process defined by the Langevin 
equation in the path integral form 

t 



dxij) 



cxp -- / dr[x + f[x[r),r)y exp - / drf{x{r),r) 

(24) ■ 

The obtained Wiener path integral (24) can be generalized for solving a system of coupled 
Langevin equations. In the next section we shall describe such a generalization. 



B. Wiener path integral for a system of coupled Langevin equations 

As a prototypical example for a system of coupled Langevin equations, let us consider a 
system of Brownian particles in random media which can be treated by the generalization 
of the Wiener path integral method depicted in the preceding section. For the system of N 
Brownian particles, there are N corresponding coupled Langevin equations 



Xi{^)+ Mx{t),t) = $i(T), 



t = l,2,...,N, 



(25) 



and equivalent ly, in the matrix form: 



x(T)+f(x(T),T)=*(T). 

As mentioned before, the next step is the functional change of variables 

y{r) = x(t) + / f (x(s), s) ds, < r < t, 
Jo 

leading to 



(26) 



(27) 



y(r) =*(t). (28) 
Similarly to the case of (19), the Jacobian of transformation can be calculated by the discrete- 



time approximation. 



J(£) = 



A{e) 

* A(2£) ... 



... * A{N6) 



(29) 
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where 



A{ne) 



1 + 



19/1(X£) 



2 dxi 

1%(X£) 

2 9a:i 



ne 



ia/i(x,£) 



2 9a:2 

1 I 1 a/2(x,£) 

"^2 9a;2 



2 axisi 



ne 



1 a/jv(x,£) 



1 + 



1 9/jv(x,£) 



2 a^i 2 dx2 ' 2 9a:;v J 

and the stars "*" in eq. (29) denote the matrix blocks which do not contribute to the 
determinant. The Jacobian of the transformation (27) is given by 



(30) 



J = hm J(e) — exp 



1^ 



2 ^ 

•4 = 1 



ds 



'0 ^-^i 

Now we can write the joint probabihty density in terms of the Wiener path integral as: 



(31) 



W^(xt,i|xo,0) 



/ 



n 



n 



dxjssij) 



X exp 



X exp 



-lE / ^T(i:, + /,(x(r),r))^ 
1 ^ 



*^^ a/.(x(r),r) 



(32) 



At this level, we have developed the Wiener path integral approach for solving the system 
of Langevin equations and we can apply this method for the case of aggregation equations. 
However, there is some conceptual point concerning the meaning of the variables in the 
aggregation equations, namely that they are not the coordinates, but are dimensionless 
velocity differences, as mentioned before. So the suggested model for the aggregation of 
particles includes the inertia of particles, but after linearizing the equations of motion one 
can get the first order differential equations in the velocity space. 

Now, it is straightforward to apply the above method to the case of two coupled Langevin 
equations. For the Langevin equations (13) describing the aggregation of particles, we can 
write the exact solution for the probability density in terms of path integral as follow: 



W^(xt,t|xo,0) 



dxAr) 



n 



dX2{T) 



X exp 
X exp 



\/47r (ir ^J^ \/ Ai^dr 

-\f^dT{{x, + hf+{x, + hf) 



dxi dx2 



(33) 



where 



fi = xi- e{Txl - xl), 
/2 = X2 - ^2 ^ + 2exiX2. 



(34) 



The obtained joint probabihty density represents an exact solution for the system of two 
coupled Langevin equations for the aggregation of inertial particles. Thus, we can write 
the maximal Lyapunov exponent in the form of a path integral, using the eqs. (15) and 



It can be seen from eq. (34) that /2 and its derivative in the expressions for W and Ai contain 
singularities at X2 — 0. However, the region X2 ~ does not contribute to the probability 
density W and the Lyapunov exponent Ai and overall the path integral is convergent. 

The case of negative Ai leads to the aggregation phenomenon. Having presented an exact 
solution of two coupled Langevin equations in terms of the Wiener path integral as in eq. 
(33), one can study all the features and properties such as Lyapunov exponent with the help 
of approximation methods, both theoretical and numerical, e.g. perturbation expansion 
in small parameters, by now very well developed lattice calculations and the saddle-point 
approximation. In this way, we can also determine the points of aggregation and non- 
aggregation phase transition and their nature, by investigating the exact expression used to 
determine the Lyapunov exponent given by eq. (35). By studying the whole integrand in the 
exponent in eq. (35) with eq. (33), as a function of xi and X2, we can find for which values 
of Xl, negative or positive, the path integral iy(xt, i|xo, 0) is larger - the region in which 
negative xi values dominate (as a function of e and F) gives the region where aggregation 
occurs. 

We will follow in this line to find analytic solutions for Lyapunov exponent by using some 
approximation methods. 

C. Evaluation of Wiener path integral for probability density in one-dimensional 
case: Instanton approximation 

As an example, we can apply the above results for the analysis of aggregation phenomena 
in one-dimensional fiuid. In the following we summarize the formulation and the results of 



(33)-(34): 




(35) 
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a specific model which has been introduced by M. Wilkinson and B. Mehlig [12'] for the 
one-dimensional fluid in a similar way as described for the three-dimensional fluid in section 
II. 

Equations of motion for any independent particle with position x{t) and momenta p{t) 

are 

x=—, p= -'yp + f{x,t), (36) 
m 



where 7 characterizes the strength of the viscous damping 12|. The random force f{x,t) is 
translational invariant in both space and time and the statistical properties of the force is 
given by the following expressions, 



=0, 



(/(x,t)/(x',t')) = «'exp 



exp 



'2^ 



(37) 



where a denotes the magnitude of the force and ^ and r are correlation length and correlation 
time, respectively. The linearized version of the equation of motion can be derived as 



6x 



6p 
m ' 



Sp = -'jSp + dn,f{x, t)6x, 



(38) 



where 6x and 6p are small separations of positions and momentums between pairs of tra- 
jectories. One can write the equations of motion in terms of X = 6p/6x (where X has a 
stationary distribution as t — ?■ 00) as the following: 



Sx = 

X = -7X - 



X6x 



m 
X^ 



m 



+ d^f{x,t). 



(39) 
(40) 



Since the separation 6x of two nearby trajectories has a lognormal distribution 
{\n\6x{t)/6x{0)\) = At, the distribution of X remains stationary. From eq. (39) we ob- 
tain for the Lyapunov exponent A 



A 



m 



(41) 



We start with the Wiener path integral analysis of this case from single Langevin equation 
(40). The probability density of particles in one- dimensional fluid, which is described by 
single Langevin equation (40), can be written in terms of a Wiener path integral. By 
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inserting / = 7X + X^/m into the eq. (24) and after restoring the diffusion coefficient D, 
we obtain 

" ± • X^ ^ ^" 

exp ^ / dr { X ^ 'yX ^ — 

X exp 



W{Xt,t\0,0) 



1 

AD 



\ 

* / 2X\ 

ctr 7 + 



m 



m J 



(42) 



with the last factor exp \ J^dr (7 + ^) 
simphfication, (42) takes the form 



being the Jacobian of transformation. After some 



W{Xt,t\0,0) = exp 



AD 6Dm 



exp 



1 



-AD I dr{X^ + ^V + 



27X3 



A dXjr) 
c(o,0;Xt,t)^J-Q VAnDdr 

x^ ADxy 



+ 



m 



m 



■)\ 



exp 



jt 
2 



(43) 



Let us consider the expression for the action in the exponent in eq. (43). 

We should mention that the normalization of the probability density is guaranteed by the 
Jacobian of transformation included in (43), although when using approximative methods 
we ought to normalize the probability function at each stage of approximation. 

The next step is to derive an analytical expression for this probability density which is 
expressed in terms of Wiener path integral (43). The eq. (43) in the present form can not be 
solved exactly, so we try to use approximative methods in some cases. Let us first introduce 
new parameters and make the following change of variables: 



f = -T = -fT,d = Dr/m% X 



-X 



-X. 



(44) 



The action in terms of new variables takes the form 



Sim - -l/lX(f)] ^ 4 /V (i(f )^ + i.t^ + A-3 + Ix' - . (46) 

where e = d/u^ = D/m?'-)^. The action in this form is suitable for WKB approximation in 
the limit e — 0. 

Here, we briefly present the WKB approximation method for the evaluation of a path 
integral (a complete description of WKB approximation can be found, e.g. in j25|). In the 
WKB approximation, a general path integral can be written by using the Taylor expansion 
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of the action 5'[X(f)] around the classical solution Xc{t), which can be obtained from the 
corresponding Euler- Lagrange equation for 



5Xr 



0, 



(46) 



where the subindex c refers to the classical solution. Using the Taylor expansion of (45) 
around the classical solution (46), gives: 



DX{f) exp 



—S[X{f) 



^exp 



—S[Xc{f) 



X 



jDr]{f)exp -^J J 



2s J J 5Xc{ti)5Xc{t2) 



(47) 



where the measures of integration and fluctuations 77(f) around the classical solution Xc{f) 
have been defined as: 



(48) 



X{f)^X,{f) + rj{f). (49) 

In eq. (47), beside the contribution of classical configurations, there are additional path 
integrals over the fluctuations 77(f) around the classical configurations. In the first approxi- 
mation, this path integral is a Gaussian one and can be calculated exactly by discretization, 
i.e. 



Drj{T) exp 



^£ J J SXc{ti)SXc{t2) 



det 



5X^ 



x=x^ 



(50) 



This path integral is only a function of time (since 77(0) — rjlu) — 0) which can be included 
in the overall normalization factor. 

Let us return to the original path integral expression (43) and use the WKB approxima- 
tion for its evaluation. As discussed above, in WKB approximation method one has to solve 

the Eulcr-Lagrange equation of motion to find the classical solutions. The classical equation 
of motion from (46) with the action in eq. (45) is : 



X,-X^- 3X^ - 2X1 + 2£ = 0, 



(51) 
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with dot denoting the derivative with respect to r. 

Instantons were introduced in the context of the turbulence for the first time by V. 
Gurarie and A. Migdal [26], through the studies of the WKB approximation for the velocity 
distribution function of forced Burgers equation. Similarly, in our problem the path integral 
can provide a framework to study the instanton contribution to the aggregation phenomena, 
although the physical meaning of the instanton solution in this case is not yet fully clear. 

Now, let us introduce a shifted variable Z = X + ^ so that the action in eq. (45) changes 



its form and the new action and Lagrangian become 



S'[Z] 







'^drdz' + liZ'^-^)' 



2eZ], 



L = ^Z^ + hz'-h^-26Z. 
2 2^ a' 



Accordingly, the Euler-Lagrange equation for general e is 



(52) 
(53) 



Z + -Z -2Z^ + 2e = 0. 
2 



(54) 



Eq. (54) in the limit e — )■ has the following instanton solution |27l |: 



Z = - tanh f - (f — fn' 
2 V2^ 



The Lagrangian for this instanton solution can be calculated and the result is 

2 



(55) 



-'inst 



1 

16 



1 - tanh^ ^ ~ 



(56) 



which after insertion in (52) gives the action 

S' = ^tanh f^(^-^o) j + ^ cosh"^ (^'^^~^°)) ^^"^^ (^*^^~^°' 



(57) 



The action (57) can be written in terms of Z,y [Z^ denotes the value of Z at the time t = v) 

as 

S' = \z, + \zj, = \z, - \zl (58) 
6 o 4 d 

using the fact that Z = ^ — for (55). By restoring the variable X and the original variable 
X, we can rewrite the action (45) as 



(59) 
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or equivalently, 



The contribution of this instanton solution to the probabihty density in the WKB ap- 
proximation (47) with (43) cancels the exponential in front of the path integral (43) and 
leads to a zero value for the Lyapunov exponent A defined in eq. (41). 

Besides the contribution of instanton solution, also anti-instanton solution which is an- 
other solution of eq. (54) can contribute to the evaluation of the path integral. The anti- 
instanton solution has the following form: 



is: 



Z = --tanh(^-(f-f^)J. (61) 
The action can be obtained in a similar way as for the instanton solution and the result 

Thus, using the eqs. (47) and (50) the WKB approximation for the probability density 
(43) yields 



ly(Xi) =/3exp 



(63) 



2D 3Dm 

From the definition of e, it can be seen that in the limit e — )■ 0, where the above solution is 
obtained, the dimensionless viscosity u is much larger than dimensionless diffusion coefficient 
d. In this regime we can consider the Gaussian part as the main part and keeping the leading 
terms in the expansion of the the remaining part, we obtain 



2D 

where /3 is a normalization constant. Hence the Lyapunov exponent can be obtained as 



(64) 



A 



1 /■+°° D 

- / dXt W{Xt)Xt = —. (65) 



We see that in the limit e — )■ 0, the Lyapunov exponent A is negative so the aggregation of 
particles (coalescence of paths) occurs in this regime. 

The obtained result A = — D/(m^7^) in (65) is identical to the one previously obtained 



m 



12| using a different approximation method. 



Along this line, one can search for corrections to the existing results by expansion around 
the instanton solution. Also the consideration of the multi-instantons configurations will be 
of interest. 
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In the following we will continue our analysis for the evaluation of the Lyapunov exponent 
in the non- aggregation phase and its derivation from the path integral formalism. 

Aggregation and non-aggregation phases can occur in many different cases for many dif- 
ferent values of diffusion and viscosity coefficients, but we restrict ourselves to special cases, 
where the path integral is exactly solvable or is appropriate for approximative methods. 
Thus, let us consider the case of a very large D, large 7 with 7/-D small but 7^/-D fixed, 
when a non-aggregation phase can be imagined, since the effect of diffusive force is more 
important than the effect of viscous force. Hence, in this case the third and fourth terms , 
i.e. /{Dm) and X'^/{Dm?) in the action in (43), can be neglected compared with the 

other terms. This brings the path integral expression (43) for the probability density into 
an exactly solvable form: 



iy(Xt,t|0,0) =exp 



7X, 



AD 



dX( 



c(o,0;Xt,t)j^ VAnDdT 



exp 



ADxy 



m 



X exp 



7t 
2 



(66) 



The expression (66) includes a Gaussian path integral, which has an exact analytical so- 



lution [25|]. Here we describe the steps of calculations very briefly. Completing the quadratic 
form in the exponent by the change of variable X — 2D/(m7^) = Y , one obtains 



M/(Xt,t|0,0) = exp 



AD 



dY{T) 

c{o,0;Yt,t)^J'Q V AirDdr 



n 



exp 



X exp 



D 



m^7^ 



exp 



2 



(67) 



After evaluation of the Gaussian path integral and restoring the original variable X, we 
reach an analytic expression for the probability density 



W{Xt,t\0,0) 



7 



AttD sinh (7t) 



exp 
X exp 



-^{X,-^)'coth{jt) 
AD ni'y^ 



AD 



exp 



D 



exp 



2 



(68) 



After taking the limit t — > 00 and normalizing the probability density W{Xt,t\0,0) in (68) 
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to unity, we obtain 



W^(Xi,t|0,0) 



7 



2txD 



exp 



D 



(69) 



Using (69) for the probability density we find the Lyapunov exponent in this case: 



A 



TO 

m 



D 



> 0. 



(70) 



Therefore, a positive quantity for Lyapunov exponent in this case for the non- aggregation 
phase is obtained. 

As a final remark, we should mention that the probability density W{Xt,t\0,0) in (42) 
satisfies the Fokker-Planck equation 

dWiXt,t\0,0) _d[{^ + 7X)WiXt,t\0,0)] d'wiXt,t\0,0) 



dt dX dX^ 

The Fokker-Planck eq. (71) can be derived by writing down the expression for infinitesimally 
shifted W{Xt+e,t + e|0,0) as in (42) and making expansion around Xt and t. Eq. (71) is 
identical to the Fokker-Planck equation used in 12|. 
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